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Abstract In this paper, we show the that the ground state energy of the one- 
dimensional Discrete Random Schrodinger Operator with Bernoulli Potential 
On , is controlled asymptotically as the system size N goes to infinity by the 

random variable £n, the length the longest consecutive sequence of sites on 
' I the lattice with potential equal to zero. Specifically, we will show that for 

O I almost every realization of the potential the ground state energy behaves 

II 1 2 . 

^ , asymptotically as (^^"^1)2 sense that the ratio of the quantities goes 

> ■ to one. 
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1 Introduction 



Random Schrodinger operators are a vast area of research in physics and 
mathematics. Since the seminal work of Anderson [2 and Mott 9 , there 
■ has been an uninterrupted fiow of ideas and publications about their spec- 

' tral properties, localization, delocalization, and dynamical properties of wave 

packets in a random potential etc. The recent experimental realization of An- 
derson localization in cold atom systems added to the interest and motivated 
new research on the nonlinear random Schrodinger equation [S] [TU] . A good 
general reference on the physics of quantum propagation in random media 
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is [T]. [6] is an excellent survey of the mathematical theory. A particularly 
important aspect of the theory is eigenvalue distribution. We refer the reader 
to [5] for a major result in this direction. 

In this paper, we focus specifically on the properties of the ground state 
of the one-dimensional Schrodinger operator with a potential generated by a 
sequence of independent, identically distributed Bernoulli random variables, 
taking values and b > 0. Studying the system in finite volume, we obtain 
detailed asymptotics of both the distribution of the smallest eigenvalue and 
the properties of the ground state eigenfunction as the system's size goes to 
infinity. Other authors have studied different aspects of Schrodinger opera- 
tors with Bernoulli random potentials. In particular, [3] contains a proof of 
exponential localization in such systems. We stress that from the point of 
view of many standard techniques, the Bernoulli case is technically difficult, 
because the distribution of the random potentials lacks density. In contrast, 
for our direct approach, the discreteness of the potential's distribution is a 
simplifying feature. It also makes it easier to compare the theoretical results 
to simulations in which using Bernoulli random variables is an advantage. 
Roughly speaking, we prove that the ground state energy is approximately 
equal to the sine wave on the longest interval on which the potential equals 
zero and control the corrections by explicit bounds. We emphasize that our 
method of estimating the ground state energy (especially from below, which 
is more difhcult) is different from that applied e.g. in [TT]. Simon's estimate 
is more generally applicable, but also less accurate while our upper and lower 
bounds coincide up to lower order terms. Accurate Lifschitz tail asymptotics, 
based on estimates similar to ones used here, are the content of a work in 
preparation with V. Borovyk. In future work with J. Stasiiiska we will ad- 
dress the analog of this result for higher-dimensional systems. [2] and [3] 
are related references on diffusions processes in potentials given by Poisson 
clouds. The work is motivated in part by research of the authors on the non- 
linear case |12j . In the remainder of the introduction, we define precisely 
the model and state our main result. Its proof relies on two bounds, which 
are proven in the following two sections. The proof is completed in Section 
4, which also includes a limit theorem for a quantity closely related to the 
ground state energy. Finally, Section 5 contains a summary and a short dis- 
cussion of further research. 

The goal of this paper is to study the ground state energy Eq of the 
operator H = —A + V on the space £^{1, N}. The Laplacian operator — A 
is considered with Dirichlet boundary conditions at the sites and N + 1 
and defined by —Af{i) = 2f{i) — f{i -I- 1) — /(« — 1) and V is an TV-vector 
with components V{j). We assume that V{j) arc i.i.d. random variables on 
a probability space {f2,J-,P) with the Bernoulli distribution: 



The ground state energy Eq is defined as the minimum of (</>, Hcj)) over 
(j) € £^{1,...,A^} with II0II = 1 (equivalently, one could take the minimum 




with probability p 

h with probability q ~ I — p 
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over (p G (^{0, 1, N,N+1} with \\(f)\\ ~ 1 satisfying the Dirichlet boundary 
conditions). We wih show that Eq is approximately the energy of the dis- 
crete sine wave S^ij) = sin( j^^^) supported on the longest island of zero 
potential. This is motivated physicaUy by the idea that finite potential bar- 
riers act essentially as infinite barriers for states with energy that converges 
to zero. 

The precise definitions are as follows. The longest island of zero potential 
is the longest set of consecutive points j € {1, ■ . • ,-/V} such that V{j) — 0. 
We denote its length £jv. There may be more than one island of length ^jv, 
but the point is inconsequential to this paper. The operator H = —A + V 
is self-adjoint, hence for every system size N and for every realization of 
the random potential there exists a ground state vector realizing the 
ground state energy. Any w £ f2 gives rise to a realization of the potential 
which defines islands where the potential variables equal zero. We denote the 
(consecutive) islands by li and the length of li by Li . The number of islands 
is a random variable denoted n (we emphasize that while A'^ is deterministic, 
n is random — this wiU be important in Section 4). It is natural to expect 
that — j> oo as — > oo with probability 1 — this is shown in Section 4. 
We can choose 'Pn to be real and nonnegative since it is the ground state 
of a nonnegative operator. (Here is the weh-known argument: the energy 
can be written as {^n.H^n) = V{j)\<pN{j)\^ + \^n{j + 1) - fjv(j)|'. If 
^Af(j + 1) ^N{j) are off by a phase factor, the energy can be made smaller 
by a complex phase factor correction. This contradicts the energy-minizing 
property of tf^Ar, therefore we can assume that 'PnU) are all of the same 
complex phase factor and we wiU take them to be real and nonnegative.) 

As a remark, our results naturally extend to the top of the spectrum: the 
highest eigenvalue behaves as b— (^i,\iyi where the roles of p and 1 —p = q 
are interchanged, so is the longest island of b potential. To prove this, we 
conjugate H with the unitary operator U : — )■ (— l)l-'l(/)(j), which trans- 
forms the Laplacian spectrum to 4 + Z\ and leaves the potential unchanged, 
so the bottom of the spectrum oi b — V ~ U~^{~A)U is the top of the spec- 
trum of —A + V. The claim now follows from our ground state result with p 
replaced by 1 — p. [11] 



Our main result, indicated informaUy above, is: 
Theorem 1: With probabihty one 

lim — = 1. 



N^oo ( 



2 Upper Bound on 

Lemma 1: For any realization of the potential V and for any system size A^, 



E^< 



'° - i^N + 1)2 
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Proof: By definition, — min{< H(j),(f) >: \\(f>\\ — 1}. To find an upper 
bound, we will construct a test function supported on the longest island, 
whose energy equals sin^ ( jj^^q-j- ) , which is of course bounded by the right- 
hand side of the claimed inequality. Restriction to functions supported on 
the longest island can be understood as a comparison to the system with 
the infinite potential barrier on the island's boundaries — see the discussion 
in the Conclusion. In the remainder of the proof, we will show that with this 
restriction the kinetic energy is minimized by SnU) = -j; sin(i^pj-). 

Consider a wave function on {0, . . . ,L} (treated as Z/(i + 1)Z), a cyclic 
group of order L + 1, see [13]) with Dirichlet boundary conditions: (/)(0) = 
4>{L + 1) = (identifying the point L + 1 with 0). Its Finite Fourier Transform 
is 

1 



VL+1 



The Finite Fourier Transformation diagonalizes the Discrete Laplacian. Tak- 
ing the Finite Fourier Transform of the Discrete Laplacian acting on a vector 
(/), we obtain 

fA^){p) = -yl=Sf^^mj) - 0(j - 1) - 0(j + 
VL + 1 



Using the fact that (/)(0) = (/)(L + 1) = and <I>(-1) = <j){L): 

1 r 2tv ipj T 2iTiv{L-^l) r 2TTip(j- 



^=L=[2 -e^- e-^]rjI-V(,)e^ 



This implies that the eigenvalues of — Z\ are sm^{^^); p = 0, . . . L. The 
smallest eigenvalue of the Laplacian without any boundary conditions comes 
from p = 0, but the corresponding eigenfunction, which is a constant, van- 
ishes because of the Dirichlet boundary conditions. Therefore, the minimal 
energy is given by p = 1 and p = L. Taking a linear combination of the vectors 

e^+i and e , corresponding to p = 1 and p = L, we obtain a nonnegative 
vector that satisfies the Dirichlet boundary conditions: 0(j) = sin(i^^). We 
calculate its norm: 



sin(^^)f = y |sin2(^^)|2 



)^|sm(-^)| ^ 
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_ L 
~ 2 

Thus the normahzed eigenvector is S'l+i(j) = •^/^sin(-^^). (Note that 

in the above we considered vectors defined on the cychc group Z/(L + 1)Z 
rather than on {0, L + 1}. However, since the discrete sine wave is zero at 
the boundary points, the inner product will not include the boundary terms 
of the Laplacian acting on the vector.) 

Wc define the test function on {0, 1, ...,7V + 1} as follows. On the 
longest island of zero potential, let = -S^iv+iO) and let = for all 
other j. The energy of this test function is 

= sin2(:^^) 



< 



JV + 1 

2 



{(n + 1)2 

The ground state energy is defined as the minimizer of < Hcj), (j) > for all 
norm one vectors. Therefore, the ground state energy must be less than the 
energy of the test function. 

n 



3 Lower Bound on Eq 

For each realization of potential and for every A'', we have a specific non- 
negative ground state In this section, we prove a lower boimd on the 
energy of \I/n, which will be our main technical tool. This bound approaches 
the upper bound of the previous section when £n ^ oo as N ^ oo. We will 
show in section 5 that this is true with probability one. We define 5]^ and 
as the values of on the sites with potential equal to b, adjacent to 
the left and the right of the island li and define rrii as the norm of !Pjv on 
li, = ll'f'jvl/JP- We will use the definition of !Z'jv as the energy minimizer 
to estimate its energy from below. As the critical tool, we first estimate its 
norm on sites of b potential, using the the upper bound on the ground state 
energy proven in Section 3. 

Corollary 1: Let 'Pn be the ground state vector for a given realization 
of potential. Let B = {j : V{j) = b}. Then ||if'7v|s|P < by Lemma 

1, where Pn\b is the restriction of S'at to the set B. 

Proof: The contribution of energy of ^'jv on B = {j : V{j) = b} is bounded 
by the quantity in Lemma 1. 

^2 
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> {V^n\b,^n\b} 

= b\\^N\Bf 

Therefore, WnIbW^ < ^ □ 

It follows from the corollary that < b(/j+iy where the 

extra factor of 2 is introduced to cover the case in which the outer boundaries 
are shared by islands. The norm of if'Ar is thus bounded below on sites of zero 
potential by \\>I^n\b4^ > 1 - bief+iy^ • 

Next, wc! fix 7 e (0, 1). While in what follows the specific value of the 
parameter 7 does not play a crucial role, it will enter the convergence rate in 
our main result and its choice may become important in further work based 
on similar ideas. 

Definition: A heavy island of is an island /j such that H^Z'jvl /< |P > 
max((5}^, 6\^)^{Li + l){£N + iy~'^ . The island which is not heavy is (naturally) 
called light. 

For the lower bound on the energy, we will estimate the contributions of 
the heavy islands of \1/n and neglect the contributions of the light islands. 
First, we show that norm of is concentrated on the heavy islands. 

Lemma 2: Let M be the index set of all heavy islands Then 

37r2 



Proof: First, we split into two orthogonal vectors, where B = {j : 
V{j) = b}. 

1 = Wn\\^ = WnIbW + Wn\b4^ 

Next, we split !i'jv|B<= into heavy islands and light islands. 

\\Wn\b4^ ^ S^eM\\^N\uf + S,eM4^Nkf 

By the definition of light islands, their contribution to the norm of !Z'Ar is 
bounded above. 

SieM4Mh\\'' < SieMcmaxidi,5}ifiLi + l){eN + l)^-^ 

< Si^M^ max((5i,(5}j)2(^jv + 1)'"^ 
< 



where Ei^Mc max(5[,^^)2 < 2||^'jv|sf < ^(ifrrTF- Thus 

1 = W^nIb^ + s.^MWMhf + ^ieM4Muf 
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Subtracting appropriate terms, the bound is shown. □ 

The lemma guarantees we have at least one heavy island when < 
1. The following theorem is the main technical result that bounds the energy 
below by estimating the contribution of energy of on the heavy islands. 

2 - 

Lemma 3: li In > ^j-'' — 1 , then the lower bound on the energy 
contributed by heavy islands is 

V 2 

[vj—^ + Oi 




and thus Eq is also bounded below by this quantity. 



Proof: We will first find an explicit form of the state iPn minimizing the 
energy on a heavy island /j (which exists by the condition of the theorem) 
using the Lagrange multiplier method with fixed nonnegative values of 5]^ 
and S]^ as well as rrii the norm of the restriction of the state to the z-th 
island. Afterwards, we will obtain information about the values of nii, S}^ 
and 6}^ in the actual ground state. 

To find the energy-minimizing vector under these conditions, consider, for 
a fixed realization of the V, the function E = {.f,Hf) = S^^^{-Af{j))f{j) 

subject to the constraint S^lif{jY = mf and /(O) = 6\ and /(£+ 1) = (5]^, 
the boundary conditions set by . Because the components of the minimiz- 
ing vector are nonnegative, we only need to consider the vectors with real 
compoments. 

Introducing a Lagrange multiplier A, the minimizing vector satisfies 

which gives the set of equations. 

dE 



dm 

dE 



4/(1) - 2(51 - 2/(2) = 2A/(1) 



dE 



4/(j) - 2/(j - 1) - 2/(j + 1) = 2A/(j) 
= Af{L)-25\,-2f{L-l) = 2\f{L) 



dfiL) 

With some rearranging, we obtain the eigenvalue equation 



-Af{j) = A/(j) 



8 



The critical points are eigenvectors of —A with boundary conditions /(O) = 
and /(-L + 1) = The set of eigenvectors contains /o(j) = + Sl 

and /fe(i) = Cfe sin(.Sfc j^^^ j + ife), where fc = 1, .... L is the frequency index, 
Ck is a normaHzing constant, and Sk and make the sine wave satisfy the 
/(O) = (5J^ and /(L + 1) = 5^ boundary conditions. Since 7, is a heavy is- 
land, llif-jvl/.f > max(J[,(5}j)2Li(£jv + 1)^"^. The square of the norm of the 
linear function /o satisfies ||/o||^ < max((5}^, and is strictly smaller 

than max((5];^, (5^)^ij(^jv + 1)^"''' , so the minimizing vector must be one of 
the sine waves. Since is strictly positive, Skk must be sufficiently close 
to 1 (we will prove a precise estimate later) and it is notationally convenient 
to put = 1, so that /i is the minimizing vector. Evaluating = 

2 2 4 4 

{-Aci sin(si;^j+ti),ci sin(si;^j + fi)) = +0( (/;+i)4 ) gives 

us the energy contribution of 'I/n \ ■ 

To estimate the energy of the ground state, we must understand how 
the coefiicients si, ti and ci depend on the boundary conditions 6]^ and 
(5^ and the mass rrii. In particular, we will find a precise estimate on si. 
The frequency index 1 is dropped from /i, si, ti, and Ci for simplicity. The 
boundary conditions are: 

csin(f) = 5}^ 
c sin(s7r + t) = 5}^ 

2 

Note that both 5}^ and 6}^ are positive and smaller than jp^^q^fp-. Since 
li is a heavy island, / cannot be monotone on all of it because in this case its 
mass there would be bounded above by max(5|^, which would contra- 

dict the definition of a heavy island. Since the ground state is nonnegative, it 
follows that s is bounded above by 1 . To solve for s, we use the left boundary 
condition to solve for t where t G [0,7r/2). 

t = arcsmf— ) 
c 

We deal with the right boundary condition similarly. As argued above, / is 

not monotone on li, so solving the right boundary condition for sir+t requires 
taking the inverse of y = sin(a;) on [7r/2, 37r/2] which is a; = — arcsin(y) + tt 

sn + t = — arcsmi — ) + tt 
c 

Subtracting the two equations we get the following equation for s. 

1 (5* 

s = 1 (arcsin(— ) + arcsin( — )) (1) 

TT c c 

To bound s from below, ^ and ^ need to be boimded above. To do this, 
calculate that the square norm of the vector restricted to the island before 
normalization. 

II sin(.s ^ a; + = Sf^^ sin^(s ^ , x + t) 
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2 



2 ^■=^ + 1 



(2) 



The second term is a Riemann sum of the function cos(2(,s7r.x + t)) muhi- 
phed by ^^^^ . The associated integral is simpler to calculate. The difference 
between the sum and integral is bounded by 



cos(2(^^^^j+t)) 



cos(2(s7ra; + t))dx 



Li + l 



(3) 



where the second term is the missing term from the Riemann sum (the right 
endpoint). We can estimate the difference between the Riemann sum and the 
integral by noting that the derivative of cos(2(s7ra; + 1)) is bounded by 2s7r. 
Therefore the expression (3) is bounded above the following error term. 



E, 



257r + 1 - 2id}if _ 2s7r + 1 - 2{d}if 



The Riemann sum is bounded below by the integral minus the error term. 
Continuing with the calculation of the square norm, expression (2) is bounded 
above by: 



Jo 



Li J 
Lj + l _ Lj + l 

2 ASTT 

U + 1 L, + l 



cos(2(s7rx + t))dx + s-k + 1/2 - {S^^f 
(sin(2(s7r + t)) - sin(2t)) + stt - (5^^ 
(2 sin(s7r + cos(s7r + — 2 sin(f) cos(t)) 



2 4s7r 

+ S7r-(^lj)2 (4) 

Now, we can substitute using the boundary equations. Since stt + f e [f , tt], 

ccos(s7r + t) = — \/c 



Therefore, the expression (4) is equal to 
Li 



1 



((5jj)2 and since t G [0, f ], ccos(t) = ^/c^ - {Si^. 

5\ 



i ^2 



-f\n- 



This is a bound for ||csin(;^^ + By definition, 

+t)f, so 



2 = ||csin(s^ 



< 



c2 + . 



STT 



(5) 



The last equation can be solved for the first term in the square brackets 
(?. We then divide by the larger of the values 5\ and Without loss of 
generality, assume (5^ < 




{Li + lW^f 



{Li + lWnf 
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Move the last term to the left side to get the inequality 

"2 / 2{sn-{6^^n 




(6) 



Since the island /i is heavy, ||/f = ||!fAr|/,f > inax(6l,6}i)^{Li + l){eN + 
and 




+ 1)1-^ (7) 

As stated before, s is bounded above by 1 and below by and the ^-terms 

are bounded above by ^J^^, so ^"""'^y ^ > -5(^7+1^^)2 > -1 which 

2 

means this term does not effect the limiting behavior of ,f, ,3 as £]v — > 00. 

Since this is a heavy island, s is strictly greater than zero because if'jv achieves 
a unique maximum on /j that is not an endpoint. The limiting behavior of 
the negative term as ^jv — ^ 00 depends only on the behavior of the square 
root terms. 

We want to show the term 2{ij^ + dominates as In — ^ 00. The order 
of ^ 




is at most the order of -n- • It follows that the order of ,^1 does not depend 
on the square root terms, so the term 2{£N + iy~^ must dominate. Therefore, 



and 



. Because (n > 1, 24"^ - ^^/^) > 4"^ and we assumed 5i < 6}^, we 

have 

(M< J_ 

C2 - ^2 -^1-7 

Since arcsin(x) < we can apply the bound to equation (1) 

1 (5* 

s = 1 (arcsin(— ) + arcsin(— )) 

TT c c 
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> 1- 

> 1 - 



c c 

2 
1 



For all heavy islands, s > 1 
each of the heavy islands is 



<l(l-7)/2 



. Then the energy contribution of 



TO,- 



[Li + I)- 



{Li + ly 



which is bounded below by 



rrii 



<,(l-7)/2 



{In + If 



+ o{ 



{In + ly 



\ 



\ 

) 

/ 



where we used < In and dropped Sj from the fourth order term since 
s < 1. The total mass of heavy islands is bounded below by ( 1 



fc(£jv + l)-' 

by Lemma 2. Thus the energy contribution of heavy islands is bounded below 

by 



377^ 

Ml 



1 



1 



«(l-7)/2 



{In + ly 



0{ 



{i 



N 



1 



Since the ground state energy is bounded below by the energy contribution 
of heavy islands, we have the desired lower bound. 

□ 

Using Lemma 1 and 3, the ground state energy is controlled as 



1 - 



1 



i^-^^'^j \{^N + l) 



{^N + ly 



) ) < < < 



{^N + ly 



Therefore, as 



•JV 



00, i?o^ 



4 Limiting Distribution of ^iv 



This section is a brief discussion of the limiting distribution of ^jv with most 
details left to the reader. It is natural to conjecture that the distribution of 
£n determines the distribution of Eq , but the bounds in the previous section 
are not strong enough to give that result. The bounds do, however, show the 
asymptotic relation 



lim 

JV->-oc 



N 



^ 1 



1 



7+1? 
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claimed in Theorem 1. We finish here the proof of this statement and, for 
completeness, we derive a limit theorem for Im- 

The first step is to find the distribution of Mj , the maximum length of j 
independent islands. This uses the theorem from [7] that states 



J 



where Xi is the individual random variable and Mj is the maximum of j i.i.d. 
random variables with the same distribution as Xi . Note that the distribution 
of an individual island is geometric: P[L = k] = p^~^q, for k = 1,2,... {k 
starts from one because an island has at least one site of zero potential). 
Because P[L < x] — P[L < [xj], where [a:J is the largest integer smaller than 
X, for the Uj{T)to have a limit we must choose a subsequence jk such that 
the fractional part of — log^ jfc, denoted {— \ogpjk}, converges monotonically 
from above or from below to 9 £ [0, 1]. When {— log^ j^} I 9 & [0, 1) we get 
the following limiting distribution: 

< u,,{r)] ^ exphp--{-+^}] + 0(1/ Jk) 

This shows the possible subsequential limits parameterized by 9. 

The second step is to study the number of islands n in a system of size A^. 
In fact, converges to = in probability. To show this, we estimate 



the probabilities of the events 



and 



by 



and 



P 



{io:n< li- -S)N\} 



{uj:n> [{- + 5)N\} 



P 



p 



{CO 



E 



+ m' > N} 



rrij + m < N} 



respectively, where mj is the sum of the lengths of one island of zero potential 
and one island of b potential and m is a correction random variable for an 
extra or missing island. We apply the exponential Chebyshev inequality to 
the latter probabilities and use moment generating functions to derive an 
exponential bound. Convergence with probability one can also be shown, 
but it is unnecessary for finding the limiting distribution. 

The final step is to estimate above and below by M^^i^^sNk and 
A/^jVfc-SATfc respectively in any realizations such that — ^| < S. The proba- 
bility that n does not satisfy this condition is bounded by an arbitrarily small 
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e provided Nk is large enough, due to the convergence in probabihty. The M- 
variables then converge (over a subsequence) to one of the hmiting distribu- 
tions found above, as A^^ — >■ oo. Note that \ogp{nNk + SNk) « \ogp{iJ,Nk) + S' , 
where 6' depends directly on 5 so the convergence of the irrational rotation 
by the {\ogp{iiNk)} is preserved. Taking S ^ wc find that the limiting 
distribution at all points of continuity (all r such that {r} ^ 1 — 0) is given 
by: 



log 



logpl 



< r 



exp 



-P 



,T-{T+e} 



This is a limit theorem satisfied by (subsequences of) £n. As a corollary, 



lim In = oo 

.N->-fx 



= 1 



which combined with our bounds on Eq gives 



lim 

JV— >-oo 



+w 

thus finishing the proof of our main result — Theorem 1. 



Remark: Attempting to derive a limiting distribution for E^'' , a change 

2 

of variables to find the distribution of ^ , ,-,2 would result in: 



log^ [nNk\ ( 7r2 TT^log^p 



27r2 log^ p 
The difference between and 



/2 



> T 



is of the order log^ [nNk\ . The cor- 



iog2 ypiNk\ 



rection term on the lower bound is at best of the order log ^ ( \_jjiNk\ ) which 
diverges when multiplied by the log^[/xA^)fcJ factor. If we try the change of 
variables with the lower bound, we will run into the same problem. This 
shows that our estimates are not sufficient to establish a limit theorem for 
the ground state energy. 



5 Conclusion 

For the one-dimensional discrete random Schrodinger operator with Bernoulli 
Potential, our key result is that is controlled in the limit as — > oo by 
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Since P [limAr_>.oo ^at = oo] = 1, it follows that 



P 




The key technical part of the proof is the derivation of the lower bound which 
uses classification of islands into heavy and light. 

As pointed out in the proof of Lemma 1, the argument used there essen- 
tially compares the given system to a system with infinite potential barriers 
rather than barriers of height b. That we are able to prove an asymptotically 
equivalent lower bound means physically that, as far as the ground state en- 
ergy is concerned, the two systems are not very different and consequently, 
its leading order behavior does not depend on b as made explicit by (8). 

It can be seen from the proof, that randomness of the potential is used 
only to show that £n — >■ oo. The argument for both the upper and the 
lower bound works for individual potential realizations which satisfy this 
condition. The same result thus holds for any deterministic potential with this 
property. One could even consider a sequence of volume-dependent potentials 
(taking values and b). For example, the theorem still holds for a sequence 
of periodic potentials whose periods diverge with N. We owe this remark to 
an anonymous referee. 

Physical intution, as well as numerical results, lead us to believe that the 
ground state is approximately the sine wave on the longest island of zero 
potential. If the ground state does localize to the largest island (or one of 
them), then the ground state energy calculation simplifies greatly since, in 
the notation of section 4, ||/|| 1 — o(l) can be substituted and it can be 
shown that s « 1 — We plan to address this issue in a forthcoming 

paper. 

We conjecture that, similar to the ground state, the low energy excited 
states will be approximated by half-sine waves and higher harmonics on the 
long islands. This may be lead to a strong form of a Lifschitz tail estimate in 
the one-dimensional discrete random Schrodinger operators with Bernoulli 
potential. Related ideas are used in work on many particle systems |12j . 
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